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In this letter I show that the recently proposed local version of the fluctuation dissipation relations 
follows from the general principle of stochastic stability in a way that is very similar to the usual 
proof of the fluctuation dissipation theorem for intensive quantities. Similar arguments can be used 
to prove that all sites in an aging experiment stay at the same effective temperature at the same 
time. 



The fluctuation dissipation relations in off-equilibrium 
dynamics are a crucial tool to explore the landscape of 
a disordered system. These fluctuation dissipation rela- 
tions are different from the predictions of the fluctuation 
dissipation theorem at equilibrium, so that sometimes 
one speaks about violations of the fluctuation dissipa- 
tion theorem. The existence of these new relations has 
been firstly proved in mean field theories Jl|-§]|, however 
it has been later shown that they follow from the general 
principle of stochastic stability (|-||] , that roughly speak- 
ing asserts that a random perturbation acts in a smooth 
way. 

The fluctuation dissipation relations can be expressed 
in a rather simple form that can be easily interpreted 
from the theoretical point of view using the concept of 
the effective temperature |^,|l0| . Moreover the main pa- 
rameters entering into the fluctuation dissipation rela- 
tions have a simple interpretation from the point of view 
of equilibrium statistical mechanics |j]J^J^] ■ This fact im- 
plies that for a given system the form of the fluctuation 
dissipation relations is universal in off-equilibrium dy- 
namics, i.e. it does not depend on the way in which the 
system is put in an off-equilibrium situation (as soon the 
system remains slightly out of equilibrium). These fluc- 
tuation dissipation relations have been amply observed in 
numerical simulations jll] jl4j and quite recently in real 
experiments 
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The case that has been mostly studied is where we con- 
sider observables that are the average over the whole sam- 
ple. In this case there is a static-dynamic relation that 
connects the dynamic fluctuation dissipation relations to 
the static average of global quantities. Recent there have 
been a few investigations on fluctuation dissipation rela- 
tions that involve only given local variables fl^,[l{| . The 
aim of this letter is to give a theoretical foundation to 
these local fluctuation dissipation relations and to de- 
rive the appropriate static-dynamic relation. In order to 
reach this goal we have to introduce a new concept: the 
probability distribution of the single spin overlap: this is 
a generalization of the well know probability distribution 
of the total overlap of the system. 

Let us recall the usual equilibrium fluctuation theo- 



rem. If we consider a pair of conjugated variables (e.g. 
the magnetic field and the magnetization) the response 
function and the spontaneous fluctuations are deeply re- 
lated. Indeed, if R eq (t) is the integrated response (i.e. 
the variation of the magnetization at time t when we 
add a a magnetic field from time on) and C eq (t) is the 
correlation among the magnetization at time zero and at 
time t, we have that R eq (t) = (3{C eq (0) — C eq (t)), where 
(3 = (kT)^ 1 . If we we eliminate the time and we plot 
parametrically R eq as function of C eq , we have that 
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The previous relation can be considered as the definition 
of the temperature and it is a consequence of the zeroth 
law of the thermodynamics. 

In an aging system the generalized fluctuation dissi- 
pation relations can be formulated as follows. Let us 
suppose that the system is carried from high tempera- 
ture to low temperature at time and it is in an aging 
regime. We can define a response function R(t w ,t) as 
the variation of the magnetization at time t when we 
add a a magnetic field from time t w on. In a similar way 
C (t w , t) is the correlation among the magnetization at 
time t w and at time t. We can define a function i? ttu (C): 
we must plot R(t w ,t) versus C(t w ,t) by eliminating the 
time t in the region t > t w , where the response function 
is different from zero. 

The fluctuation dissipation relations state that for 
large t w the function R tw (C) converge to a limiting func- 
tion R(C). We can define 



dR 
dC 
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where X(C) = 1, for C > = lim^oo C eq (t), and 
X{C) < 1 for C < Coo. 

The shape of the function X(C) gives important infor- 
mation on the free energy landscape of the problem, as 
discussed at lengthy in the literature |l|-^|. It has been 
shown that in stochastically stable systems the function 
X{C) is related to basic equilibrium properties of the 
system. Let us illustrate this point by considering for 
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definitiveness the case of a spin glass. Given two equilib- 
rium configurations a and t, we define the global overlap 
as 



=l,iV ' 



AT 
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where iV is the total number of spins. The function Pj(q) 
is the probability distribution of the overlap for a given 
sample and the function P(q) is defined as the average of 
Pj(q) over the samples. 

It is convenient to introduce the function x{q) defined 
by dP(q)/dq = x(q). The relation among the dynamic 
fluctuation dissipation relations and the statics quantities 
is simple X(C) =x{C). 

There are recent numerical results jl8],[l9) that indicate 
that the fluctuation dissipation relations and the static- 
dynamics connection can be generalized to local quanti- 
ties in systems where a quenched disorder is present and 
aging is heterogeneous pp[ . 

For one given sample we can consider the local inte- 
grated response function Ri(t w ,t), that is the variation 
of the magnetization at the point i at time t when we add 
a magnetic field at the point i starting at the time t w . 
In a similar way the local correlation function (7, (t w , t) is 
defined to be the correlation among the spin at the point 
i at different times (i.e. t w and t). Quite often in system 
with quenched disorder aging is very heterogenous: the 
function Cj and Ri change dramatically from one point 
to an other |po| . 

It has been observed in simulations (lljjnj that local 
fluctuation dissipation relations seems to hold 



dR % 
dC~ 



(4) 



where Xi(C) has quite strong variations with the site i. 

It has also been suggested that in spite of these strong 
heterogeneities, if we define the effective /3f at time t 
at the site i as 



dRi(t w ,t) 
dCi(t w ,t) 



= f3Xi{t w ,t)=pf f (t w ,t) , 
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the quantity (t w , t) does not depend on the site [jl9| . 
In other words a thermometer coupled to a given site 
would measure (at a given time) the same effective tem- 
perature independently on the site: different sites are 
thermometrically indistinguishable. 

These empirical results call for a theoretical explana- 
tion. The aim of this note is to show that these results 
are consequence of stochastic stability in an appropriate 
contest and that there is a local relation among static 
and dynamics. The crucial step consists in defining a 
local probability distribution of the overlap for a given 
system at point i (i.e. Pi(q)). We will then introduce 
the notation Xi(q) = Pi(q')dq' and we will show that 
the static-dynamic connection for local variables is very 



similar to the one for global variables and it is given by 
X i (C)=x i (C). 

It is evident that we need a definition of the local over- 
lap that should be rather different from the usual one. 
Indeed the local overlap of two equilibrium configura- 
tions (i.e. (TiTi) is always equal to ±1. If we use a naive 
definition, we find that the probability distribution of the 
local overlap is the sum of two delta functions at ±1. 

An useful definition of the local overlap an of its proba- 
bility distribution can be obtained if we consider M iden- 
tical copies (or clones) of our sample: we introduce NxM 
of variables where a = 1, M (eventually we send M to in- 
finity) and N is the (large) size of our sample (i — 1, N). 
The Hamiltonian in this Gibbs ensemble is just given by 



H K {<J)= J2 H(a a )+eH R \a] 



(6) 



a=l,M 



where H(a a ) is the Hamiltonian for a fixed choice of the 
couplings and the Hr{ct\ is a random Hamiltonian that 
couples the different copies of the system. A possible 
choice is 



H R ia] = 



E 

= l,M;t= 



(7) 



l,N 



where the variables Kf are identically distributed inde- 
pendent random Gaussian variables with zero average 
and variance 1. We can consider other ways to couple the 

systems (e.g. H R [a] = E a ,b=i,M ; i=i,N K ?' b <<r b i) but lct 
us stick to the previous one for definitiveness. 

Our central hypothesis is that all intensive self average 
quantities are smooth function of e for small e. This hy- 
pothesis is a kind of generalization of stochastic stability. 
According to this hypothesis, when e — > 0, the local cor- 
relation function and the response functions tend to the 
result for e = uniformly in time. 

In the case of small, but non-zero e, let us take two 
equilibrium configurations a and t. For given K, we 
consider the site dependent overlap 



qi(a.r) 
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We define the AT-dependent probability distribution 
P t K (q) as the probability distribution of the previous 
overlap. By averaging over K at fixed e, we can define 



P?(q) = P t K (q) , 



(9) 



where the bar denotes the average over K. We finally 
define 



fta) = limi?(«) 



(10) 



where the limit e — > is done after the limits M — > oo 
and A^ — > oo (alternatively we keep e 2 M and e 2 N much 
larger than I). 
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The proof of the local fluctuation dissipation relations 
can be obtained by copying, mutatis mutandis, the proof 
of the usual fluctuation dissipation relations || . For ex- 
ample we can consider the following perturbation 



AH 



(2) _ 



Eua.b „a „b 



(11) 



a=l.M.b=l,M 



where the variables h are Gaussian random variables with 
zero average and variance S/M. If we assume for simplic- 
ity a Langevin type of time evolution, the same steps of 
give : 



X f (t)=2 / dt'Ci{t',t) 



, ^dRi{t',t) 



dt' 



The previous equation can be rewritten as 



X f\t) = 2 I dCCXi(t){C) , 



where we have defined 



dRi(t',t) fdCi(t',t) 



dt' 



dt' 



(12) 



(13) 
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If we make the crucial hypothesis that the limits t — > oo 
and 5 —> may be exchanged, we get 



lim X f\t) = XT' = /3 / dqP t {q)(l - <?) . (15) 



Generalizing the previous arguments we get 



X^(t)=s I dt'Q(t',t) 



dt' 



dCC^XiWiC) 13 / dqPi(q)(l - q s ) = (16) 



sfi J dq Xi(q)q s 1 

Using the previous equations for all values of s, we arrive 
to the conclusion that the quantity Xi(t)(C) has a limit 
when the time goes to infinity (i.e. to formulation of the 
local fluctuation dissipation relations). The limit is given 
by 

Xi(C) = lim Xi{t){C) = Xi {C) = [ dqPi(q)dq . (17) 

t-»oo J Q 

The previous equations give the connection among static 
quantities and the local fluctuation dissipation relations. 

The thcrmomctric indistinguishability of the sites can 
be proved by considering two far away sites i and k and by 
introducing a perturbation that depends both the spins 
at i and the spins at k. 

A typical example is 



AH 



(3,2) _ \p 
i,k / j 
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where the variables h are Gaussian random variables with 
zero average and variance <5/M 4 . If we proceed as before, 
we find that the probability distribution Pi,k{qi,qk) can 
be written in terms of Xi(t) and X\~{t). By imposing that 
Pi,k{qi,qk) is a positive distribution, we finds that 



Xi(t) = X k (t) 



(19) 



This conditions is just what is needed to impose the ther- 
mometric indistinguishability of the sites during aging. 
As a by product we obtain 



Pi,k{qi,qk)= I dx5(qi - qi(x))5(q k - qk(x)) (20) 
Jo 

where qt(x) is the inverse function of Xi(q). This relation 
has some interesting consequences that will be investi- 
gated in details elsewhere. 

We note that the probability distribution of the lo- 
cal overlap Pi(q), being related to a dynamical quantity, 
must depend only on the local environment around the 
point i and therefore it must have a straightforward limit 
when the volume of the system goes to infinity (e.g it 
should be independent of the boundary conditions). It is 
remarkable that for far away points (i, k) the two proba- 
bility distributions Pi(q) and Pk(q) are independent one 
from the other, but the joint probability distribution of 
qi and qk does not factorize has shown by eq. |2^. 

Summarizing, we have seen that for a given sample it 
is possible to give a definition of a local probability distri- 
bution of the overlap, that depend on the site. The prop- 
erties of this local probability distribution are related to 
the local fluctuation dissipation relations, that automat- 
ical follows from the present formalism. The condition of 
thermometric indistinguishability of the sites turns out to 
be a byproduct of our approach: during the aging regime 
all the sites are characterized by the same effective tem- 
perature. 

It is important to stress that here (as well as in the 
original paper || ) in the all the derivations we have made 
no esplicite assumption concerning the validity of static 
or dynamical ultrametricity or on the existence of a large 
separations of the time scales. Although the generalized 
fluctuation dissipation relations have been firstly derived 
in ultrametric mean field models , their validity relies 
only on more general properties like stochastic stability. 

I am grateful to Jorge Kurchan for an illuminating dis- 
cussion. 
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